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Abstract 

The optical conductivity is studied numerically for the ferromagnetic metal- 
lic state close to the charge ordering observed in perovskite manganites. We 
focus on roles of the charge fluctuation on the incoherent component in the 
optical spectra, and consider the Coulomb interaction in the one- and two- 
dimensional spinless fermion models. The numerical results show that the 
incoherent part is not noticeable, and near the quarter-filling, most of the 
spectral weights are transferred to the Drude component. In the one dimen- 
sional model, there is little or no spectral weight in the energy region between 
the Drude weight and the remnant of the absorption edge which is pushed up 
with the carrier doping, in contrast to the two dimensional one. The discrep- 
ancy between our calculations and the experiments suggests that the orbital 
degrees of freedom play important roles. 



After the discovery of High-Tc cuprates, the electronic structures in many transition- 
metal oxides have been reinvestigated intensively from the renewed viewpoints. One of 
the attractive materials is manganese oxides with the perovskite structure Ai_ x B x Mn03 
(A = La, Pr, Nd, B = Sr, Ca) owing to the negative colossal magnetoresistance (CMR) and 
its potential applications |T]-§[. The parent compound, LaMnOa, is an antiferromagnetic 
insulator where the valence of the Mn ion is +3 and the electronic configuration is 

(^2p) 3 (e g ) 1 because of the strong Hund coupling. It is noted that, as well as the spin and 
charge, the e g electron in Mn 3+ ion has the orbital degrees of freedom. The layer-type 
antiferromagnetic spin alignment and the two-dimensional character of the spin wave 
are attributed to the orbital ordering accompanied with the Jahn- Teller lattice distortion 
10-^. As the concentration of the divalent cation is increased, the system turns to the fer- 
romagnetic metallic (FM) states at x ~ 0.175 in Lai_ x Sr x Mn03, and CMR is observed near 
the ferromagnetic Curie temperature (T c ). With further doping, T c is gradually decreased 



and the charge ordered (CO) phase comes out around x ~ 0.5 in Ndi_ x Sr x Mn03 |T0[ and 
La!_ x Ca x Mn03 In the CO phase, the spin and orbital orderings coexist. Therefore, it is 
recognized that the perovskite manganites are one of the typical materials where the spin, 
charge and orbital degrees of freedom are strongly correlated. 

In the FM phase, some anomalous features in the optical conductivity cr{ui) are recently 
reported ||11|| . u{uj) shows the two-component feature, that is, the sharp Drude peak with 
the small intensity and the broad incoherent band up to about leV. Even at the lowest 
temperature, where the magnetization is almost saturated, the integrated intensity of these 
components is still growing up with decreasing the temperature. Because the spins are 
perfectly aligned at that temperature, the results indicate that another degrees of freedom 



are still active. In our previous paper ||12| , we focused on the orbital degrees of freedom 
in perovskite manganites and investigated the nature of the orbital fluctuation in the FM 
phase. With considering the strong electron-electron interactions and the anisotropy of the 
electron transfer, it is concluded that the long range orbital ordering is prevented by the two- 
dimensional character of the orbital fluctuation. We proposed the orbital liquid states where 



the orbital fluctuation is still active down to low temperatures. The numerical calculation 
based on this scenario shows the broad spectra in the optical conductivity. It has been 
proposed that the orbital liquid state is one of the promising candidates for the origin of the 
incoherent spectra observed in perovskite manganites. 

As well as the orbital degrees of freedom, the Mn ion in the FM phase has the charge 
degrees of freedom, as was pointed out above. In the CO phase around x ~ 0.5, the 
charge excitation gap, showing the gentle slope, is found in the optical conductivity in 



Pro.6Cao.4Mn03 [fL3|l . Introducing the charge carriers or applying the magnetic field to 
the CO state, the charge ordering easily melts and it turns into the FM phase at low 
temperatures. In the FM phase near the charge ordering, the charge fluctuation is expected 
to be large due to the competition between the electron kinetic energy and the long range 
Coulomb interaction, and may bring about the incoherent component. Therefore, it is 
important to study to what extent the remnant of the charge ordering is relevant to the 
optical spectra. 

In this letter, we study numerically the optical conductivity in the FM phase near the 
charge ordering. We focus our attention to the roles of the charge degrees of freedom, because 



the orbital degrees of freedom have been investigated in our previous paper JT2J . We employ 



the exact diagonalization of the finite size systems, and conclude that its contribution is less 
dominant on the incoherent broad band in the region x < 0.4 of our interest. 

We consider the e g electron with the single orbital in the perfectly spin polarized ferro- 
magnetic state. In order to describe the fermionic excitations in this phase, we introduce the 
spin- and orbital-less fermion in each Mn sites. The on-site Coulomb interactions between 
the e g electrons and the Hund coupling between the e g and t2 9 spins are irrelevant, while 
the inter-site Coulomb interaction dominates the charge ordering and its fluctuations in this 
case. We consider the following Hamiltonian, 

H = tJ2 (cfa + H.c) + > (!) 

<ij> <ij> 

where cj is the spin- and orbital-less fermion operator at site i, and rij is the charge density 



operator defined as rii = c\ci. t and V is the transfer integral for the fermion and the 
Coulomb interaction between the nearest neighbor sites, respectively. In order to obtain 
the ground state wave function and the optical conductivity spectra, we utilize the exact 
diagonalization method i.e., the Lanczos method and the recursion method in the finite 



size system |jl4] , l5|] . With considering the symmetry of the lattice and the capacity of the 
memory in the computer, the one dimensional lattice with Nl = 16 x 1 sites and the two 
dimensional one with Nl = 4x4 sites are studied. The size effects are studied in the one 
dimensional lattice, and the convergence of the results is confirmed in the 16 x 1 sites case. 
The periodic boundary condition is adopted to avoid the fictitious low energy excitation 
near the free boundaries. 

In general, the optical conductivity is decomposed into the two parts, i.e., the Drude 
component and the regular part, as follows, 

a{u) = D5(cu) + a re9 (to) , (2) 

where a re9 (u) is defined by 

a M = 7rL ^r 6 ^ ~ E n + E ) . (3) 

n ^ &n ~ ^0 

In Fig. 1, we present the doping dependence of the regular part of the optical conductivity 
spectra in the one and two dimensional cases. The Drude part will be discussed later. In the 
case where the ground state is degenerate, the spectra are averaged in the momentum space. 
The value of V/t is chosen as 12 and 4 in the one and two dimensional case, respectively, 
in order to true up the absorption edges. Although these values are taken to be larger in 
comparison with that in the actual compounds, this is for the purpose to avoid the finite 
size effects. The fermion concentration n = N/Nl is identified with x in manganites, and 
n = 0.5 corresponds to the half-filling. In Fig. 1, the optical gap is clearly shown at 
uj ~ V and at uj ~ 3V in the one and two dimensional case, respectively. This is nothing 
but the charge ordered insulator, where the charge alternation occurs in the one and two 
dimensional lattices. The optical gap corresponds to the charge transfer gap, because an 



transfer of the fermion requires the energy of V and 3V in the one and two dimensional 
cases, respectively. With decreasing the number of the fermion, the system departs from 
the commensurate charge ordered phase. From the results of the calculated Drude weight, 
it is concluded that the system is metallic, except for the half-filling case. As the system 
departs from n = 0.5, in the one dimensional case (Fig. 1(a)), the energy at the absorption 
edge is gradually increased with the reduction of its intensity. There is little or no spectral 
weight between the Drude part and the remnant of the absorption edge above V, and the 
spectral weight is transferred from the later to the former with doping. This picture for 
the optical spectra in the one dimensional chain are consistent with results predicted by the 



bosonization procedure ||16|| , where the above two components in the spectra is interpreted 
as the intra- and inter-band transitions at the Luther- Emery point. Near the quarter-filling 
(n = 0.25) in Fig. 1(a), the Drude part governs the almost all spectral weight and there are 
no remarkable weight in the finite energy region. On the other side, in the two dimensional 
4x4 lattice, the charge gap is collapsed and the spectral weight appears within the gap 
with doping, in contrast to the one dimensional case where the remnant of the absorption 
edge is pushed up. This is because the charge excitations at u> ~ V and u ~ 2V become 
possible by doping. However, we cannot find any remarkable spectra in the region of u < V 
in Fig. 1(b). At the quarter- filling, only the small amount of the spectra remains in the 
finite energy region similarly to the one dimensional case. 

Next we study the doping dependence of the Drude weight. In Fig. 2, we show the 
Drude weight (D) and the total spectral weight (/), in one- and two-dimensional cases. The 
Drude weights in the finite size cluster with the periodic boundary condition are numerically 



estimated by the following relation [17,18 



e 2 ~ N L e*N L ^ E n - E ' 1 ) 

where the first term in the r.h.s is the total spectral weight /, and < T > is the expectation 
value of the kinetic energy in the ground state. The validity of the above formula has been 



carefully examined in the one dimensional Hubbard model by the several authors fll7| , |18 



It has been confirmed that for iV — > oo, D estimated by Eq. (^) is close to the Bethe- 
ansatz value, and its doping dependence is qualitative correct even in the small cluster 
system. At half-filling in Figs. 2(a) and (b), we find small positive and negative values of 



D being attributed to the finite size effects. As the several authors have mentioned |]IT,IB 
it becomes inconspicuous with increasing the size and the magnitudes of the interaction. 
With doping into the insulating states at the half-filling, the Drude weight appears and the 
system becomes metallic. In the one dimensional case, D is rapidly growing up with the 
carrier density 1 0.5 — n|, and its value approaches to /, immediately. With further doping, 
both D and / go toward zero at n — 0.0 and 1.0 cases. On the other hand, in the two 
dimensional lattice, D is suppressed near the half-filling, reflecting the spectral weight in 
the finite frequency region, as shown in Fig. 1(b). As the system approaches to the quarter- 
filling case, however, almost all the weights go to the Drude component. With increasing V/t, 
the rising of D and / around the half-filling become slower, and the Drude part dominates 
the total spectra more and more. 

In Fig. 3, we present the static charge correlation function Xc{k) calculated in the 
finite size system. Xc(k) is defined by Xc(k) = F.T. < SriiSrij > with 5rii = rii — n. At 
half-filling, Xc{k) has a sharp peak at it /a and (7r/a, 7r/a) in the one and two dimensional 
lattice, respectively, reflecting the commensurate charge ordering. With doping carriers, the 
intensity of XcC 71 "/ ) in the one dimensional lattice is drastically decreased and the broad 
incommensurate peak appears around k = ir/a — 5. In the cases of n > 0.4375 (more than 
one hole doping), the magnitude of Xc(k) is smoothly reduced with the doping. On the other 
hand, in the two dimensional 4x4 lattice, the noticeable peak at (jr/a,ir/a) still remains 
up to the case of n = 0.375. With further doping, the intensity of the correlation function 
is almost saturated, and the incommensurate peak appears at (7r/2a,7r/a) at the quarter 
filled case. The doping dependencies of Xci^/a) and xd' K / a i it /a) are strongly correlated 
with these of the Drude weight, shown in Fig. 2, both in the one and two dimensional cases, 
i.e., as the charge correlation at the boundary in the Brilloun zone is diminished, the Drude 
weight is gradually growing up. It is suggested that the strong charge correlation prevents 
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the coherent motion of the doped carrier. 

As we show in Figs. 1 and 2, any conspicuous spectral weight are not found in the region 
of to < V in both one and two dimensional cases. This feature is quite in contrast to the 



results for the positive-U Hubbard and t- J models pl9|-|2T|j which has been extensively studied 
in connection with the mid-infrared absorption observed in High-Tc cuprates. As the carriers 
are doped in the half-filled antiferromagnetic insulator, the spectral weight appears around 
uj ~ J ~ t 2 /U, in addition to the u ~ U structure originated from the transition between 
the Hubbard bands in the Hubbard model. In the two dimensional case, this structure 
becomes remarkable and its weight is comparable to that located around U. The origin 
of this mid-infrared component is interpreted as the incoherent motion of the doped holes 
which creates distortions in the antiferromagnetic background. In a word, the presence of 
the spin degrees of freedom is essential to the emergence of the band. On the contrary, 
only the charge degrees of freedom is taken into account in the present model, describing 
the perfectly polarized ferromagnetic states without the orbital degeneracy. There are not 
additional degrees of freedom which govern the optical spectra smaller than V. 

Turning to the experimental results, the frequency dependence of the conductivity spec- 
tra is measured in Lai_ x Sr x Mn03 at x=0.0 (antiferromagntic insulator) |22|| , 0.1(ferromag- 
netic insulator), 0.175(FM) and 0.3(FM) |ll|]. The weight and width of the incoherent band 
becomes the most intense at x=0.175 and it entirely fills up the optical gap (about leV) ob- 
served at x=0.0 and 0.1. With further hole doping (x=0.3), the weight is rather diminished 
in spite that the charge fluctuation is expected to be large, as the system approaches to 
x=0.5. The incoherent band also appears by applying the magnetic field in Pro.6Cao.4Mn03 



where the charge ordering is realized without the magnetic field. ||13|| . At about 6.5T, the 



spectra are continuously spread up to about 0.8eV and entirely covers up the charge gap 
observed at zero magnetic field. On the contrary, in the present theoretical calculation, 
the almost all weights are transferred to the Drude part near n = 0.25 where the charge 
correlation at the zone boundary are fairly weaken, and even around the half-filling, there is 
no noticeable spectra in u < V. Although the numerical calculations are limited in the one 



and two dimensional finite size systems, we have no doubt that the contributions from the 
charge fluctuation are less dominant on the incoherent component observed in the perovskite 
manganites. 

The results presented in this paper suggest the other candidates for the origin of the broad 
incoherent component in the optical spectra, such as the low dimensional orbital fluctuation 



proposed by the present authors ||12|| , the interband transition proposed by Shiba. et.al. 



|1 and so on. Actually, incoherent band based on the orbital fluctuation are calculated by 



the static approximation |12] where the orbital fluctuation is treated as the slowly varying 



degrees of freedom. It is shown that the width of the spectra is an order of the transfer 
integral and its weight is expected to be growing up with decreasing the temperature. Also, 
our previous calculation is consistent with the observed two-component feature in the optical 
conductivity with considering the interaction between the orbital and Jahn- Teller distortion. 
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Figure captions 

Fig. 1: The optical conductivity spectra in the one dimensional 16 sites (Fig. (a)) and 
the two dimensional 4x4 sites (Fig. (b)). V/t is chosen as V/t = 12 and V/t — 4 in the 
one and two dimensional case, respectively. The calculations are performed at u + i5 with 
S/t = 0.0002. 

Fig. 2: The Drude weight (D/e 2 ) and the total spectral weights (ir < T > /N) in the one 
dimensional 16 sites (Fig. 2(a)) and the two dimensional 4x4 sites (Fig. 2(b)). 

Fig. 3: The charge correlation function Xc(k) = F.T. < 5ni5rij > in the one dimensional 16 
sites (Fig. 3(a)) and the two dimensional 4x4 sites (Fig. 3(b)). The parameter is chosen 
as V/t = 4 in the both cases. 
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